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ABSTRACT 


Richardson's "four-thirds law" of horizontal diffusion 
was tested using aerial photography as a data gathering tech- 
nique. Plywood floats and current crosses suspended both near 
the surface and at nine feet were used as diffusers. The 
scales investigated ranged from 10 to 525 meters. The inves- 
tigation was conducted in 36 fathoms of water, 3000 meters 
from the nearest land in Monterey Bay, California. Stommel's 
(1949) method of analysis was used. 

The results indicate”a clear” dependence of diffusion on 
diffuser weight and lend some evidence to Robert's (1961) 
theory of turbulent diffusion, in that the diffusion increases 
more rapidly with scale than proposed by Richardson (1926). 
This conclusion is supported by the use of confidence limits 


upon the data. 
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1. Horizontal Turbulent Diffusion and Richardson's 


"Four-Thirds Law" 


Oceanic diffusion is not well understood and a diffi- 
cult process to describe mathematically. The one-dimen- 


Sional Fickian law of diffusion is 


ot a x 
where Vv is the concentration of a substance, ZX is posi- 
tion, t€ is time and K is the diffusivity (a measure of the 
rate of diffusion) of the substance. This equation ade- 
quately describes molecular diffusion, which is dominated 
by concentration gradients. Using the Fickian equation, 
Stommel (1949) pointed out that the probability of two 
discrete particles initially a distance De apart being a 


distance by apart at a time t later would be 


(b,-b.)> 
Plb.b) = saeee exe b ERY 





In effect this equation states that the probability of the 
diffusing particles' change in separation with time is a 
maMetLoOn of (b = b,)* only, and not DA Gx by: This result 
is not supported by experimental data. 

Richardson (1926), on the basis of six atmospheric 
diffusion experiments having length scales of from 1.5 x 


ie. cm. EOn > <x 10° 


cm., postulated a new equation to describe 
diffusion in a turbulent medium. He noted that the basic 
difference between molecular diffusion and turbulent diffu- 


sion was that the independent variable for the latter was 


the particle's separation from its neighbors and not the 


particle's position, and proposed the equation 
O te 2 = S) 


Here, { is the distance between the particles and is called 
the “neighbor separation." The number of particles having 
neighbors with neighbor separations between L and gee 2 
1s given by ¢(Lld 2 where q ( L) is the "neighbor concentra- 
tion." The neighbor concentration can also be regarded as 
the probability that two particles will be a distance L 
apart. The term calet is the "neighbor diffusivity" and 
is anale gous to the diffusivity in the Fickian equation. 
Atmospheric experimental data and mathematical simplic- 


ity led Richardson to suggest that F( 2) has the form 


Pid) KAY? 
where W is an empirically derived constant. This is 
commonly referred to as the "four-thirds law." Stommel 
(1949), on the basis of data obtained using aerial photoqg- 
raphy, suggested that the four-thirds law was also applica- 
ble to oceanic turbulent diffusion. If it is assumed that 
the physical governing quantities are L and @ (the tur- 
bulent energy dissipation per unit mass) dimensional analy- 


4% 1/3 
sis shows that F = c 2 E 


, where C is dimensionless. 
Thus K= ce'l2 , so that K will vary with the turbulent 
energy dissipation, or probably the turbulent energy itself. 


This is certainly reasonable and leads one to expect a 


variety of values of K 
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Dis The Standard Method of Determining Neighbor Diffusivity 


as a Function of Neighbor Separation 


The commonly accepted method of determining oceanic 
neighborhood diffusivity as a function of neighbor separation 
is by means of pairs of floats released simultaneously on the 
water surface. Let Le be the separation of a vair of floats 
at time ty and J the separation at a later time ty: The san 
terval of time ti - ty is chosen so that the ratio ( D- fs)/Le 
is of the order of 1/10. The time interval is then a function 
of the scale, Ly . Richardson's equation may then be written 


as 








a 
which has the solution 


Se val an iz 


Thus the a sens has a Gaussian distribution. 


st = 


The standard deviation of ie from the ees is f2tF( 9.) 


so that the neighbor diffusivity for the value Lo is 


Pays oe 
ries 





where the bar denotes the mean of many pairs. Since all pairs 
will not have the same initial separation BU , a mean value of 
g is used: 7) = (&- f)” 
L . ce 
Many individual studies have been conducted to test 
Richardson's four-thirds law. These are listed in the refer- 
ences. While general agreement exists as to the validity of 
the four-thirds law, the data has never been sufficient to 


actually confirm the law, and a value of the constant C has 


i bal 


not been determined. Rather, a range of values of K can be 
found in the literature? trome6-00G. towed, og Bed ae, These 
values have usually been found by passing a best-fit line of 
4/3 slope through the data points plotted on a log-log scale. 
Little attention has been paid to the slope itself, or to 
random sampling errors and resulting confidence eae on both 
the slope and the constant K. This should be corrected, 
especially as recent developments in turbulence theory 


(Roberts, 1961) predict that the slope should be 3/2 instead 


of 4/3. This field study is a first step in that direction. 


she Objectives of the Field Study 


The field study has two main objectives. First, to in- 
vestigate Richardson's four-thirds law using as many observa- 
tions as practical. Second, to further investigate sub- 
surface horizontal diffusion as was first done by Snyder 
(1967), who found a marked change in diffusivity with depth. 
The greatest possible scale commensurate with technique and 
available equipment was studied. Measurements of scales 


ranging from 10 to 525 meters were obtained. 


4, Techniques Employed to Obtain Diffusion Data 


A T-ll Fairchild aircraft mapping camera was obtained 
from the Naval Oceanographic Office, Suitland, Maryland. This 
camera uses a 6-inch £:6.3 (Bausch and Lomb) class T metroqon 
lens to obtain 9-1/2 x 10-1/4 inch negatives from a film sup- 
ply roll with a 390 foot capacity. This lengen of seein awe 


sufficient for between 450 and 460 negatives. The operation 


VZ 


of the camera is completely automatic and is driven by a 27.5 
volt d.c. motor which derives its power from the aircraft's 
main electrical system. A vacuum of 2" mercury is used to 
hold the film flush against the focal-plane frame. The inner 
edges of this focal-plane frame form a 9 x 9 inch perimeter 
for a picture area on each film negative. Other features in- 
corporated in the T-ll camera (which make it ideally suited 
for diffusion studies) are an altimeter, a navigational clock, 
geclaeoy recorder tand a built-in filter wath a 1.5 _andex.. ~Time, 
altitude and other pertinent information such as date of 
flight were thus recorded on the margin of each negative. 

The sweep second hand allowed time to be read to within 0.01 
seconds. Comparison of the aircraft and the camera altimeters 
showed them to agree perfectly. The camera was mounted in a 
Navy US2-A aircraft in such a manner that level flight en- 
sured the axis of the camera to be in the vertical. From the 
aircraft's flight instruments it was ascertained that all 
picture-taking runs were within two degrees of level flight. 
Operating under these conditions ensured a linear length 
scale over each negative for easy measurement. 

The diffusing particles consisted of 2' x 4' x 1/4" ply- 
wood floats, weighing about 10 pounds. To each float was 
attached a metal current cross consisting of four mutually 
perpendicular vanes, 6 inches wide and 35 inches long. The 
current crosses weighed nine pounds, and required additional 
buoyanc*’” to remain afloat. Styrofoam sheets were used for 
this purpose and an almost neutral buoyancy was achieved. 


This minimized any direct effect of the wind on the floats. 
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In the first three runs half the floats had current crosses 

at nine feet; the remaining floats had current crosses just 
beneath the surface. All floats were painted white. A 10- 
inch red stripe was painted diagonally across the floats whose 
current crosses were near the surface. These red stripes were 
easily recognizable on the negatives. A total of twenty-three 
floats was used. 

Data were taken on March 21, 28, April 4 and April 18, 
1968. The final run made used seventy-four floats without 
current crosses. A 40-foot motor launch set out the diffusers 
on the first two runs and a 63~foot converted torpedo re- 
triever on the last two runs. 

The floats were put out in pairs; the individual floats 
composing each pair were 30 feet apart. On the first three 
runs each pair consisted of a float with a current cross at 
nine feet and one with a current cross near the surface. The 
initial pattern was laid out in a circle with an approximate 
diameter of 300 feet. The diameter of the initial pattern on 
the fourth run was 450 feet due to the large number of floats 
used. After laying the pattern the boat stayed as near to the 
center of the pattern as was possible without causing propel- 
ler wash (which could have a definite harmful effect on the 
study). 

The experiment was conducted in Monterey Bay, California, 
in approximately 36 fathoms of water and about 3,000 meters 
from the nearest point of land. The exact initial position 


each time was 36°38.6' north latitude, 123°53.3' west longitude. 
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The pictures were taken at intervals of between three and four 
minutes at altitudes of 1,000 feet for the first three runs 
and 2,000 feet on the fourth run. Wind conditions were taken 
as those reported by Monterey Municipal Airport, and are shown 


in Table I below. 


TABLE I 
Wind Data 
Run Race Wind Direction Wind Velocity 
iH 3-21-68 100° - 140° 55.0 m/sec 
2 3-28-68 250° -= 280° 2.5 - 4.0 m/sec 
3 4-04-68 320° - 340° 1.5 - 3.1 m/sec 
4 4-18-68 300° = 320° 2.5 - 5.0 m/sec 


During the second run a swell of four to six feet from 
the west (the same direction as the wind) was reported by the 
motor launch. No significant wave or swell activity was en- 
countered on runs one and four. Some breaking waves were no- 
ticed from the air during the early stages of run three, but 
these subsided before the run was completed. The sighting of 
whitecaps and the report from the motor launch of actual sea 
conditions was not commensurate with the reported wind, which 


is therefore suspected to be in error. 


oe Data Analysis 
A series of one to three pictures at three-second inter- 
vals was taken on each pass of the aircraft over the diffusing 


particles. Without the use of optical pnhoto-navigational 


Bs) 


equipment, and using only the motor launch as a reference sys- 
tem, it was very difficult to fly directly over the center of 
the diffusing particles. As a consequence some floats did not 
appear in all pictures and were therefore eliminated from 
further analysis. Table II contains the pertinent information 


for each of the four runs. 


TABLE II 


Run 1 (23 floats) 


Altitude 1 Altitude 2 Delta t Pass 
(feet feet Pini in i 
1,000 1,180 4.09 1&2 
1,180 e210 4.14 Za 
a ale 1,060 se 3&4 
1,060 1,160 Ar ARO 4&5 
ly 1.6.0 1,060 3.92 5&6 
1,060 14020 4.44 6&7 
dy 02,0 1,080 ore yi 7&8 
1,080 1,040 4.50 8&9 
1,040 1,000 4.36 9&10 
1,000 1,040 4.31 10&11 
1,040 1,060 4,39 11812 
1,060 Tg 0.0.0 4.40 12813 
1,000 1,060 Ses 13814 

Run 2 (14 floats) 
10 980 2153 1&2 
980 1,020 3.92 2&3 
ete 1,020 Avlige3O 6&7 


Altitude 1 


feet 
i O20 
eos 0 
Jer ON es 
i008 
980 


960 


NOTE: 


1,060 
1,080 
Ole 
1,040 
1,080 
1,080 


1,040 


NOTE: 


2,040 
2,000 


i 2.0 


Toes Li eteconmt.) 


Alerteade. 2 
(feet) 


1,080 
1,000 
1,000 
980 
960 


260 


Passes 4 & 5 were not used 


Delta t 
min 


5.6.5 
S267 
B07 
AO 
320 


Sie 


Runs (cl etoats) 


1,080 
p00 
1,040 
Hier el o319, 
i, SO 
1,040 


1,000 


4.36 
4.10 
8.54 
8.00 


A ee 3 


Run 4 (44 floats) 


2,000 
iy, Je0 


2,000 


ay 


Pass 


7&8 

8&9 

9&10 
10&11 
eo alge 


P2643 


because less than 


half the floats appeared in each of the negatives. 


1&2 
2&3 
3&4 
4&5 
7&8 
10&11 


11&12 


Passes 6 & 9 were not used because less than 


half the floats appeared in each of the negatives. 


L&2 


2&3 


3&4 


A transparency was made for each negative, and was used 
with a Travel-Graph overhead projector. An enlargement factor 
of five was obtained when the transparencies were projected 
on a plain white paper sheet stretched and taped to a wall. 

A simple color coding scheme was used to mark each float on 
the paper sheet. Each pass of a run was represented by its 
own symbol and each float was further identified as to vertical 
location of the current cross. In this manner the position of 
each float was tracked throughout the duration of the run. 

The paper sheet was placed on a cork bulletin board and a 
colored pin placed in the center of each mark representing 

the position of a float. A different color pin was used for 
each pass. This simple scheme allowed for positive identifi- 
cation of each float. The ranges and bearings of the floats 
were then measured from a convenient reference point near the 
center of the pattern. This reference point could vary from 
pass to pass. Ranges were measured to within 0.025 inches 

and bearings to within 0.25 degrees. A digital computer was 


used to calculate initial distance J. , changes in distance 


d.-@> | mean distance ~ (@+%.) , and neighborhood diffusivity 
" 
oe for each neighbor pair. Ranges measured in inches 


were converted to the proper scale in meters using the rela- 
tionship; 


IMAGE GROUND COVERED 


FOCAL LENGTH ALTITUDE : 
Histograms were drawn for each run and values of + ( Ves 2.) 
were separated into class intervals. The sizes of the class 


intervals were selected in such a manner that each group had 
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about the same number of values. The values of F in each 
class were then averaged to obtain the F for that class 


Incerval. 


oS: Results and Conclusions 
A total of 2,374 observations covering scales of 10 to 
525 meters were made. The results are given in Tables III - V. 
TABLE III 


{current crosses near surface) 


ily, 


Run 1 

@lass Interval Number of il E( 0) 
(meters) Values Avg. (meters) (meters sec +) 

Os = eG | 47] 3 170 0.0104 
Shor “79 46 70.6 020392 
Go =) 198 44 AEP 0.0628 
98 =- ll2 43 LEL..0 0.0500 
112 - 146 42 132.20 OZ 55 
146 - 170 61 153.0 O21 GLY 
f70-= 198 45 Ge 0 OLA TSS 

Run 2 
10 - 61 18 36.0 0.0096 
61 - 92 on 76.4 ORO el ves: 
92 - 107 23 9822 0.0602 
EOF = 3122 23 a0 0.0884 
122 - 174 one a2 70 0.1040 


Class Interval 


(meters) 
6 12 2 
Wee lB 3 
ines 83 
183 200 
200 228 
228 253 
Pais: 305 
305 B52 


TABLE III 


Riana 


Number of 


Values Avg. 


30 


2.5 


OAL 


25 


2a 


a7 


20 


22 


20 


(GOtii=:) 


ie 


(meters) 


65D 
legal. 0 
164.0 
O30 
23 20 
2500 
2on.; 0 


33970 


F(X) 
(meters“sec” 
0.027 
OGIO 
0.0854 
02510 
O.n2daed 
0.1490 
0.1 7aeo 


0.4140 


1 


) 


TABLE IV 


(current crosses at nine feet) 


Run 1 
Class Interval Number of DQ FQ) 
(meters) Values Avg. (meters) (meters“sec +) 

10 - 61 oe. 40-7 0.0130 
61 - 79 49 ere: O20393 
79 - 98 47 eRe rers: 0.0464 
98 - 116 48 EOS 02 O20778 
116 - 134 52 1262s 021096 
134 - 153 a7 141.9 0.1374 
SG) Seem 4g Gall 78 174.4 0.2584 

Run 2 
21 - 61 25 43.4 0.0146 
61 - 92 35 Thee 0.0510 
eee) 107 23 syed: OL One USyS, 
oy — 137 25 Ae) 0.0655 
37 — 256 29 184.0 O21923 

Run 3 
10 - 107 42 764-2 0.0356 
107 - 152 39 B29 0 0.0750 
152 - 183 48 166.0 0.0948 
183 - 228 47 20720 0.1450 
228 - 276 35 247.0 0.3460 
276 - 350 45 balay C, 0. 5060 


Zi 


TABLE V 


(no current crosses) 


Run_4 

Class Interval Number of pe Ee 
(meters) Values Avg. (meters) (meters“sec +) 

10 - 91 OE eile 0.0226 
Siva eZ 85 124.5 0. L073 
F52 = 202 I/O. Les O72945 
212) -=5200 105 254.0 0.4475 
290 - 320 107 308.0 0, Gore 
320 - 348 106 334.0 0.7760 
348 - 385 P22 364.0 0.8560 
385 - 425 ey 402.0 0.9248 
425 - 525 114 452.0 1.0780 


The results were also plotted on logarithmic paper, as 
neighbor separation L versus neighbor diffusivity F( 2) 
(Figures 1-3). The values of K varied little from run to 
run even though whitecaps were present on one of the runs 
and wave and swell conditions varied from run to run. 
Table VI gives the values of K for visually estimated 
best-fit lines of four-thirds slope and compares these 
values with those of Ozmidov (1957) and Snyder (1967). 
These values are accurate to only one significant figure 


owing mainly to sampling errors associated with confidence 


limits divseussed below. 


a2 


TABLE VI 


(values of K in a See) 

Figure 1 (current crosses near surface) W = 0.003 + 0.002 
Figure 2 (current crosses at nine feet) K = 0.003 + 0.002 
Figure 3 (plywood floats without 

current crosses) Ao = 0 0062 on 02 
Snyder (1967, paper sheets) Koe= 02006 
Snyder (1967, current crosses at nine 

feet) Ko = 02002 

Szmidoy (1957) Ko O00 ee 0 007 


NOTE: Snyder did not report confidence limits on K 


The values of K do not change appreciably with depth. The 
weight of the diffusers, however, does effect the value of 
K by a factor of about two. Snyder's (1967) result using 
current crosses was lower than any calculated in this study. 
This may be due to either sampling errors or possibly long 
term variations in € . Using diffusers without current 
crosses yielded a value for K that is in agreement with 
those of both Snyder and Ozmidov (1957). 

Eighty-percent confidence limits were calculated for 
F( 2) using the number of values averaged for each point as 
the degrees of freedom of a chi-square distribution. This 
procedure is based on the fact that Stommel's approximate 
solution of the neighbor diffusion equation shows le) 
to be normally distributed. While a line of four-thirds 
slope fits the plotted data points quite well, in all cases 


a line of steeper slope gives a better fit. For example, 


ZS 


Loled 


the relation F(L)= K 2 gives the best-fit line for 

Figure 3. Snyder, using paper sheets as diffusers, obtained 
PYVA 

a best-fit relation F(L)= KR for data obtained in the 


same area as this experiment. 

Roberts (1961) theory of turbulent diffusion, based 
on Kraichnan's (1959) work, predicts a 3/2 slope. The 
results all suggest a power law fractionally larger than 
the one studied, thus offering some experimental confir- 
mation of Robert's theory. However, this is certainly not 
conclusive, and warrants further investigation. To do this 
a study of diffusion on widely separated scales would be 


useful. 
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PLYWOOD FLOATS WITHOUT 
CURRENT CROSSES 





ie eye slope 
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1a) 
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FIGURE 3 
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APPENDIX £ 


Description of Camera Installation 

The T-1ll aerial mapping camera was installed in a 
Navy US2-A aircraft just aft of the radome. The radome and 
associated electronic equipment had been previously removed, 
so the weight of the camera (110 pounds) did not present 
a stability and therefore safety- problem. 

The access door for the radome compartment was removed 
and an aluminum plate fashioned to the same size as the 
access door set in its place. The nenae and latch arrange- 
ments of the original installation were used to hold the 
aluminum plate in place. A hole of identical size and shape 
as the lens housing was cut into the plate to accommodate the 
camera. Four, twelve-inch steel bolts were fastened perven- 
dicular to the aluminum plate with aircraft self-locking 
nuts. A steel retaining ring was placed over the steel 
bolts and held in place by leveling nuts. The pvnlane was 
then positioned in a level flight attitude by means of 
hydraulic jacks, and the camera secured to the plate by the 
retaining ring. The retaining ring leveling nuts were ad- 
usted so that the lens axis was set in the vertical. This 
critical alignment was made with the aid of the liquid level 
bubbles on the camera housing. 

Electrical power came directly off the airplane's main 
d.c. power supply. The source of the two-inch mercury 
vacuum was the de-icing system. Circuit breakers, located 


in the cockpit, were installed to isolate the camera's 
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electrical system as an added safety precaution. The 
Circuit breakers were kept open at all times except when 
actually engaged in picture taking runs. During these runs 
the anti-collision light located immediately forward of 


the camera lens was turned of£ for obvious reasons. 
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